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Abstract. 

Analysis of high-resolution stellar spectra relies heavily upon atomic data. These 
include energy levels, wavelengths, cross-sections for various types of interactions be- 
tween particles and photons, such as photoionization and collision induced transitions. 
Quantum-mechanical calculations provide largest part of these data. In this paper I de- 
scribe atomic data necessary to compute model atmospheres and line formation for the 
late-type stars both in the assumption of local thermodynamic equilibrium (LTE) and in 
a more general case of non-LTE. I will focus on transition metals with 21 < Z < 29 and 
discuss whether (and where) more complete and/or accurate atomic data are necessary. 



1. Introduction 

Chemical composition of late-type (FGK) stars provides major observational constraints 
to the chemical evolution of the Galaxy. These stars span a huge range of metallicities 
and ages and, due to relatively cool surfaces, their spectra are rich in lines of different 
atomic and molecular species that allows to determine abundances of many chemical 
elements. For certain elements, lines due to forbidden and dipole-allowed transitions, 
as well as molecular lines, can be analysed simultaneously, e.g. [CI], CH, and C I. 

At present, determination of element abundances usually proceeds in two steps. 
First, one constructs a model of a stellar atmosphere using basic conservation laws and 
certain simplifying assumptions, which are valid for a restricted range of stellar pa- 
rameters. The model gives variation of basic atmospheric quantities, local kinetic gas 
temperature T e and number density of free electrons N e , as a function of depth. Next, 
the model is used to solve radiative transfer equation for individual spectral features 
to compare them with observations. In both steps, it is necessary to describe correctly 
absoiptive properties of the atmospheric plasma, which depend on the distribution of 
atoms and ions among excitation states and ionization stages Nj, and on the cross- 
sections for interaction of gas particles with radiation field. In a simple case of local 
thermodynamic equilibrium (LTE), the former is computed from the Saha-Boltzmann 
equations for local values of T e and N e at each point of the model atmosphere. How- 
ever, in reality atomic level populations Ni are, in return, affected by the radiation field, 
which is generally highly non-local due to scattering processes. Thus, Nj also depend 
on physical conditions at other depth points and must be determined from the solution 
of full equations of statistical equilibrium for each ion and each of its energy levels, an 
approach known as non-LTE. These equations must include the rates of all processes, 
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which represent interaction between particles and photons, such as radiatively-induced 
excitation and ionization of atoms. The need for accurate cross-sections of these pro- 
cesses renders atomic data as one of the key ingredients in stellar atmosphere studies. 

Below, I will briefly describe atomic data necessary to compute model atmo- 
spheres and line formation both in LTE and non-LTE cases. I will focus on the late- 
type stars with 4000 < T e ff < 6500 and 3 < logg < 5, and on transition metals with 
21 < Z < 29. Abundances of these elements in old stars provide key information for un- 
derstanding explosive nucleosynthesis in supernova e. Brief reviews o f ato mic data used 
in pre cision stellar spectroscopy are also given in IJohanssonl (120081) and Mashonkina 



(2009). 



2. Opacity in model atmospheres and spectrum synthesis 

Construction of theoretical model atmospheres requires a detailed knowledge of the 
monochromatic absorption coefficients over a huge frequency interval. For the range 
of stellar parameters corresponding to late-type stars, the following opacity sources are 
relevant. The continuous extinction is represented by the free-free transitions in H i, 
He , H~, and H^, Rayleigh scattering for H, He, and H2, and Thomson scattering on 
electrons. Also included are bound-free transitions in H~, H 1, He 1, H~ and H+. Very 
important absorbers in the UV and visual spectral range are C 1, Mg 1, Si 1, Al 1, Fe i. For 
these species it is crucial to have accurate photoionization cross-sections. iGruppI (120041) 



showed that quantum-mechanical cross-sections for Fe I (IB autistal 1 1 9971. Iron project) 
lead to a much better agreement between theoretical and observed solar fluxes in the 
UV and IR spectral regions. The influence of Fe I bou n d-free opacity on the solar UV 



continuum fluxes was also demonstrated in iBell et al.l (120011) using MARCS and AT- 
LAS9 models, although they argued that the quantum-mechanical cross-sections should 
be increased by a factor of two to match the fluxes in the spectral range 3000 . . . 4000 
A. These same arguments may not apply to other types of atmospheric models. 

Line opacity is dominated by bound-bound transitions in neutral and singly-ionized 
atoms of Fe-peak elements, although opacity due to other atoms (e.g. H 1) and molecules 
(e.g. H2, CH, OH, Mg H, TiO) has to be in cluded as well. Modern model atmosphere 



codes (e.g. MAFAGS: iGrupp et all 120091) use tens of millions of lines to construct 



opacity distribution functions or sample them at certain frequencies to save computa- 
tional efforts. Recent calculations of KurucJ] increased the number of predicted radia- 
tive bound-bound transitions for several neutral atoms, such as Fe 1, Ti 1, Cr 1 to few 
millions; most of these tran sitions originate from the levels with very high excitation 



energies. IGrupp et al.l (120091) showed that including all predicted transitions has a small 
effect on the mean temperature structure of the solar atmosphere: up to 25 K in the 
line formation regions. Another interesting result was a reduced discrepancy between 
observed and synthetic solar fluxes in the UV. 



3. Line formation 

Calculation of synthetic profiles for comparison with observations is not trivial even 
under LTE and requires atomic data for individual lines: wavelengths, energies of 
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Figure 1. (a) The line of Co n in the solar KPNO flux spectrum (filled circles). 
The profile computed with A = 49 mK for the lower level a 5 P 3 gives an excellent fit 
of the observed profile (black trace). The profiles computed neglecting HFS (dashed 
trace) or with the less accurate laboratory value A = 40 mK (dotted trace) are also 
shown, (b) The solar line of Fe i computed with various scaling factors to inelastic 
H i and e~ collisions, Su respectively S e . The observed profile is shown with filled 
circles. 



lower levels of transitions, oscillator strengths, pressure broadening parameters, hy- 
perfine splitting and isotopic shifts. For most of these parameters, data from laboratory 
measurements are usually available, however the accuracies are not always as good as 
desirable. 

In calculations of abundances by the most reliable and rigorous method of spec- 
trum synthesis, one derives the product of abundance and oscillator strength of a tran- 
sition, gfe. Thus, in absolute abundance analyses of stars the accuracy of former is 
never better than the latter. In g eneral, it is more dif ficult to measure accurate oscillator 
strengths for weak transitions dNilsson et al.ll2003l) : this is manifested in the fact that 
a typical error in experimental /-values of weak lines is of the order of 20% or more. 
Often one has to rely on oscillator strengths from different sources, affected by various 
sources of uncertainties. 

In analyses of stars other than the Sun, there is a way to avoid ^/-values com- 
pletely. In the relative abundance analyses, oscillator strengths do not enter calcula- 
tions at all, because any abundance estimate derived from a single line in a spectrum 
of a metal-poor star is referred to that from same line in a spectrum of a reference star, 
e.g. the Sun. However, this works well on condition that for both objects a line appears 
on the same part of the curve-of-growth and, thus, is equally sensitive to various atmo- 
spheric quantities (abundance, pressure and turbulent broadening). This is true only for 
a very limited range of stellar parameters. 

Some lines are affected by hyperfine splitting (HFS) and/or isotopic shift, and these 
types of line "broadening" are comparable or even larger than the thermal broadening 
of lines at local kinetic temperatures in the atmospheres of solar-type stars. HFS is par- 
ticularly pronounced for odd-Z elements Mn and Co with large atomic masses and large 
nuclear spins, 5/2 respectively 7/2. The effect of HFS is to de-saturate spectral lines 
leading to a compound profile, where the strength of each component is linearly propor- 



4 



Maria Bergemann 



tional to the element abundance. To compute accurate separation of HFS components, 
magnetic dipole constants A must be known with an accuracy of at least 10%. Fig. Q] 
shows that the Co 11 line at 3501 A computed w ith an approximate estimate A = 49 mK 
for its lower level a 5 P 3 dPickering et al.|[l998 ) cannot fit the observed profile at all. In 
contrast, a more accurate laboratory value A = 40 mK from iBergemann et al. (2010) 
provides an excellent fit to the observed profile. A profile computed without HFS is 
also shown. 

Fin ally, accurate wavelengths of transitions are imp ortant in studies of line asym- 
metries ( Dravins et al.lll98ll : iNordlund & Dravins|[l990 ). e.g. in diagnostic analyses 
by means of 3D hydrodynamical model atmospheres and for identification of blending 
features. 



3.1. NLTE radiative transfer 

Computing NLTE line formation is very time-consuming and complex compared to 
LTE, because much more atomic data are necessary in addition to the standard set of 
parameters for individual lines as described above. 

Atomic models for statistical equilibrium calculations must be fairly complete. 
Various processes representing interaction between atoms, electrons, and radiation are 
to be included. These are photon absorption in line transitions, photoionization, ex- 
citation and ionization by collisions with free electrons and neutral hydrogen atoms. 
All processes include their reverse reactions; for scattering in discrete radiative tran- 
sitions, one usually assumes complete frequency redistribution. The current status of 
NLTE modelling for neutral atoms of transition metals is described below. These ele- 
ments have a partly filled 3d subshell that gives rise to fairly complicated but interesting 
atomic properties. 




Figure 2. Grotrian diagram of the Fe i atom constructed using the new predicted 
levels and transitions from the Kurucz's database. The ground state of the Fe n, a 6 D, 
is indicated. 
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The number of levels and discrete radiative transitions in neutral atoms of the Fe- 
group is enormous. Recent calculations of Kurucz predict tens of thousands of levels 
and few millions of transitions. Such atomic models are not tractable even with ID 
NLTE codes; thus, recently, we developed efficient algorithms to combine atomic levels 
and transitions into super-levels and super-lines in a way that a reduced model atom 
inherits the basic properties of the complete model and has similar performance under 
restriction of different stellar parameters. For example, application of these algorithms 
to Fe i reduced the number of levels from 37 500 to 292, and the number of radiative 
bound-bound transitions from - 6029000 to 13600 (Fig. |2]). Furthermore, for 3D 
NLTE calculations, we have constructed a model with 172 levels and 1120 transitions, 
which is now being tested in ID and 3D NLTE line formation codes. The requirements 
to the accuracy of transition probabilities, or /-values, are not very strict. Strong lines, 
which are subject to non-equilibrium excitation effects, have laboratory /-values with 
5 - 10% accuracy. For transitions between a bulk of intermediate- and highly-excited 
levels, the accuracy of 20-30% is sufficient, because the energy separation of the levels 
is in most cases small enough for collisions to dominate over radiative rates. 

The situation is more uncertain for bound-free radiative processes. Hydrogenic ap- 
proximation is usually ad opted for photo ionizatio n; however co mparison with quantum- 
mechanical data for Fe i ( Bautistalll997 ) and Cr i ( Naharll2009l) shows that the former is 
not correct even within an order of magnitude. The quantum-mechanical cross-sections 
for the levels of both atoms are larger in the background and are characterised by promi- 
nent resonances that leads to significantly enhanced ionization rates in statistical equi- 
librium calculations and, thus, to different abundances. In some cases, accurate position 
of resonances is important, because UV radiation field varies strongly with frequency 
in the atmospheres of late-type stars and ionization rates may change depending on 
whether a resonance appears at energies with high or low fluxes. Available laboratory 
data very sparse; e.g. for a few Ti I intermediate-excitation levels the cross-sections 
were measured by lYang et al.l ( 2009h by the method of resonance ionization mass spec- 
trometry. Experimental values are very useful to check the accuracy of theoretical val- 
ues for individual transitions, however they cannot cover the complete atomic system 
that is essential in NLTE calculations of atomic level populations for stellar atmosphere 
studies. Thus, quantum-mechanical data are indispensable. 

Cross-sections for H i and e~ collisions are basically unknown. Thus, we have 
to rely on commonly-used approximations or to neglect these processes completely. 
Cross-sections for continuum and allowed discr ete trans i tions due to inelastic colli- 
sions with H i are computed with the formulas of iDrawinl ( 19681 fl969allbh . This recipe 
was originally developed for col lisions between equal at o ms, H and Ar, and later it was 
extended to Li-H collisions by ISteenbock & Holwegerl ( 19841) . The rates of allowed 
and forb idden transitions due to col lisions with e lectrons are calculated from the for- 
mulae of Ivan Regemorteri (Il962h and lAllen1 ( 19731) . res pectively. Bou nd-free transitions 
caused by electron collisions are treated according to ISeatonl (119621) . Similar to pho- 
toionization, the accuracy of the data obtained with these f ormulae is very low. Com- 
parison of R-m atrix calculations for e~ collisions for Ca II ( Melendez et al.|[2007h with 
the formula of Ivan Regemorteri (1 1962b reveals order of magnitude differences. What 
concerns inelastic H i collision, compared to the Drawin's formulae, ab initio quantum- 
mechanical calculation s predict significantly lower collision rates for ce rtain transitions 
of simple alkali atoms ( Belyaev & Barklemll2003l ; lBarklem et al.ll2010l) . and they show 
that, in addition to excitation, other effects like ion-pair formation become important. 
Excitation and ionization balance in the atoms of the Fe-group depend on e~ collision 
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efficiency for solar-type stars, i.e. metallicities [Fe/H] > -0.5. In the atmospheres of 
metal-poor stars, where NLTE effects on abundances are large, collisional excitation is 
fully controlled by neutral hydrogen. 

The deficiencies in NLTE models often reveal themselves, when synthetic spectra 
are compared with observations, e.g. by producing large abundance discrepancies be- 
tween different spectral lines of the same atom. Thus, it is common to introduce scaling 
factors, which are adjusted to obtain agreement between various spectral lines. Fig. 
\Vp shows the profiles of the solar Fe i line computed with different scaling factors to 
Drawin's H i collision rates, S h = , 0. 1 . Inelast i c e~ c ollision rates computed using the 
formulae of Ivan Regemorterl ( 1962 ) and ISeatonl (ll962l) are also scaled by S e = 0.01, 1. 
The figure shows that a proper choice of S h and S e may also remove discrepancy be- 
tween observed and synthetic line profiles. 



4. Discussion and outlook 

The calculations of NLTE stellar spectra for the late-type stars are usually decoupled 
from calculations of model atmospheres. As a rule, the latter are computed assum- 
ing LTE, and line formation calculations are performed keeping the input model fixed. 
Such a "simplified" NLTE approach is not fully justified, because various elements af- 
fect the structure of the atmosphere by contributing free electrons and/or opacity, such 
as Mg, Al, Si, Fe, Ti (see Sec. [2]). Deviations from LTE in the excitation and ion- 
ization equilibria of these elements will lead to changes in the energy balance. Thus, 
NLTE rate equations for all these elements must be solved simultaneously in model at- 
mosphere calculations, a numerically challenging task. Another problem is that atomic 
data for some of these elements are of very low quality or lack completely. This refers 
to photoionization, H i and e~ impact excitation and ionization cross-sections. 

The first problem has been successfully overcome with enormous progress in nu- 
merical techniques and computing facilities in the pas t 10 years. At present it is possible 



to compute NLTE line-blankete d model atmospheres ( Andersonlll989l : IShort & Hauschildt 



2005al : lHaberreiter et all [2008 ). These and other studies proved the importance of 



NLT E effects on the atmospher ic structure and spectral energy distributions for the 
Sun (IShort & Hauschildt 2005a ) and non-negligible differences for metal -poor stars 



(IShort & Hauschildtll2005bh . There is no doubt that soon self-consistent NLTE models 
will replace LTE-based ones; accurate atomic data will stimulate and, thus, accelerate 
this transition. 
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